Di-electron spectrum at mid-rapidity in p + p collisions at y/s = 200 GeV 
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We report on mid-rapidity mass spectrum of di-electrons and cross sections of pseudoscalar and 
vector mesons via e + e~ decays, from ^/s = 200 GeV p+p collisions, measured by the large acceptance 
experiment STAR at RHIC. The ratio of the di-electron continuum to the combinatorial background 
is larger than 10% over the entire mass range. Simulations of di-electrons from light-meson decays 
and heavy-flavor decays (charmonium and open charm correlation) are found to describe the data. 
The extracted u — > e + e~ invariant yields are consistent with previous measurements. The mid- 
rapidity yields (dN/dy) of <f> and J/tp are extracted through their di-electron decay channels and are 
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consistent with the previous measurements of <f> — > K + K and J/ip — > e + e • Our results suggest a 
new upper limit of the branching ratio of the r\ — > e + e~ of 1.7 x 10 -5 at 90% confidence level. 



PACS numbers: 25.75.Dw 

I. INTRODUCTION 

Di-leptons are a crucial probe of the strongly inter- 
acting matter created in ultra-relativistic heavy-ion col- 
lisions [H [|| . Leptons are produced during the whole evo- 
lution of the created matter and can traverse the medium 
with minimal interactions. Different kinematics of di- 
lcpton pairs (mass and transverse momentum ranges) 
can selectively probe the properties of the formed matter 
throughout its entire evolution 0, H| . 

In the low invariant mass range of produced lepton 
pairs (Mu < 1-1 GeV/c 2 ), vector meson in-medium 
properties (mass and width of the spectral functions 
of p(770),w(782), and (/>(1020)) may be studied via di- 
lcpton decays and may exhibit modifications related to 
possible chiral symmetry restoration [H, HJ . For example, 
at the SPS, an explanation of the low-mass di-lepton en- 
hancement in the CERES e + e~ data from Pb+Pb colli- 
sions requires substantial medium effects on the p-meson 
spectral function Q. Also, NA60 recently reported a sig- 
nificant excess of low-mass /i + yT pairs in In+In collisions 
above the yield expected from neutral meson decays |6| 
which is consistent with a broadened spectral function [7[ 
but not a dropping-mass scenario @. 

At RHIC, the PHENIX experiment observed a signif- 
icant enhancement for 0.15 < M ee <0.75 GeV/c 2 in the 
low transverse momentum (px < 1 GeV/c) part of the 
di-electron continuum in Au+Au collisions compared to 
that expected from hadronic sources [9j . Models that suc- 
cessfully describe the SPS di-lepton data consistently fail 
to describe the PHENIX data in the low-mass and low-pr 
region dGil. Also, in the higher px range, direct photon 
yields were derived through di-electron measurements at 
RHIC allowing an assessment of thermal radiation [ll[ . 
Additional precision experiments with large acceptance 
and a broad range of beam energies can provide invalu- 
able insights in this subject [l|. 

The di-lepton spectra in the intermediate mass range 
(1.1 < Mu < 3.0 GeV/c 2 ) are expected to be directly 
related to the thermal radiation of the Quark-Gluon 
Plasma (QGP) @,I3. However, significant background 
contributions from other sources have to be measured 
experimentally. Such contributions include background 
pairs from correlated open heavy-flavor decays, which 
produce a pair of electrons or muons from the semilep- 
tonic decay of a pair of open charm or bottom hadrons 
(cc — > l + l~ or 66 — » l + l~). In the high-mass region 
(Mu > 3.0 GeV/c 2 ), J/tp,T, and their excited states 
are used to study the color screening features of the 
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QGP [12]. The PHENIX collaboration has reported di- 
electron spectrum in p+p collisions and has found the 
data are very well described by the hadronic cocktail 
and heavy flavor decays for the entire mass range within 
the uncertainty of the data and the cocktail 9]. The 
first di-electron continuum measurement from STAR in 
y/s = 200 GeV p + p collisions, presented in this pa- 
per, provides a crucial reference for corresponding future 
STAR measurements in heavy-ion collisions. 

Rare processes like leptonic decays of hadrons provide 
possible observables to be used in searching for traces of 
new physics Beyond the Standard Model (BSM) (T^H^ . 
These decays usually involve electromagnetic or weak 
couplings which can be calculated to a high degree of ac- 
curacy within the Standard Model (SM). In addition to a 
direct observation of the Higgs boson, the Large Hadron 
Collider (LHC) looks to explore BSM physics. Deviations 
of rare process observables from SM predictions may be 
taken as indirect evidence of a new coupling beyond the 
SM physics [H} , which can also be explored at lower ener- 
gies. The pseudoscalar mesons (for example, 77 or 77' ) are 
particularly interesting since their decay to e + e~ pairs is 
suppressed by a 2 (10 -4 ) and by helicity conservation due 
to the small electron mass (r 2 = (m e /m v ) 2 ~ 10~ 6 ) p7| . 
The branching ratio (B.R.) of ?/ — > e + e~ is 2.3 x 10~ 9 ac- 
cording to the SM predictions, however, couplings from 
BSM physics may increase this B.R. significantly [l3j . 
RHIC offers high-luminosity nucleus-nucleus collisions 
with large multiplicities and copious hadrons of inter- 
est thereby providing a unique environment for study- 
ing rare decay processes, nuclear medium effects, and 
searching for BSM physics. With recent upgrades, in- 
cluding the new Time-Of-Flight detector (TOF) [3 and 
improved data acquisition system [l9j . the STAR exper- 
iment is able to benefit from a high rate capability as 
well as excellent lepton identification at low momentum 
in the search for rare decays. 

This paper is organized as follows. Section [TT] shows 
the detector and data sample used in this analysis. Sec- 
tion IIIII describes the analysis details including electron 
identification, electron-pair distributions, background 
subtraction, and di-electron continuum without efficiency 
correction. Section ITVl presents the details of the simula- 
tions of di-electrons from light-meson decays and heavy- 
flavor decays, collectively called cocktail simulations. 
The efficiency correction for the di-electron continuum, 
the corrected di-electron continuum, and systematic un- 
certainties are also discussed in this section. Results on 
the yields of ui, </>, and J/ip from di-electronic decays are 
presented in detail in Sec. El The rare decay of 77 — > e + e~ 
is discussed in Sec. EH Lastly, Sec. I VIII provides a con- 
cluding summary 
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II. DETECTORS AND DATA SAMPLE 

Two sub-detectors are used in this analysis at mid- 
rapidity at STAR [2(| : the Time Projection Chamber 
(TPC) [U and a newly installed TOF [H|. The TPC 
is the main tracking detector at STAR, measuring mo- 
menta, charge, and energy loss of charged particles. The 
TPC, immersed in 0.5 Tesla solenoidal magnetic field, is 
a 4.2 m long cylinder surrounding the beam pipe. Its 
fiducial volume ranges from 50 to 200 cm radially and is 
filled with P10 gas (90% argon and 10% methane). Elec- 
trons from ionized gas drift toward the nearest endcap, 
where the signals are read out. The TPC readout is di- 
vided azimuthally into 24 sectors, 12 at each end. Each 
sector is divided into inner and outer subsectors with 
a total of 45 pad row readouts. The pad row readouts 
provide precise positions of primary ionization clusters 
(hits) generated by charged particles. The ionization en- 
ergy loss of a charged particle in the TPC gas (dE/dx) is 
used for particle identification [H, [23[ ■ The most proba- 
ble dE/dx is determined from the mean of a distribution 
in which the 30% of clusters with the largest signals are 
discarded (i.e. a truncated mean). For the data taken in 
2009 and analyzed here, 72% of the full TOF system was 
installed and operational. The full TOF system contains 
120 units which we call trays, 60 in the pseudo-rapidity 
range < rj < 0.9 and 60 for —0.9 < rj < 0, with each tray 
covering 6 degrees in azimuth. The TOF has a typical 
stop timing resolution of 85 ps allowing the identification 
of tt(K) up to a momentum of 1.6 GeV/c and p(p) up to 
a momentum of 3 GeV/c pLI^. 

The minimum-bias triggered events were defined by 
the coincidence of signals in the two Vertex Position De- 
tectors (VPDs) [H, located on each side of the STAR 
interaction region and covering 4.4 < \rj\ < 4.9. This 
di-electron analysis used 107 million minimally biased 
events from non-singly diffractive (NSD) Jli = 200 GeV 
p + p collisions (a NSD = 30.0 ± 3.5 mb [27|), in which 
the collision vertex is required to be within 50 cm of the 
mean of the distribution, nominally at the center of the 
TPC and along the beam line. 

III. DATA ANALYSIS 

A. Tracking and particle identification 

Hits belonging to charged particles traversing the TPC 
are collected and reconstructed into tracks with well de- 
fined geometry, momentum (p), and dE/dx. Only tracks 
that project back to within 1 cm of the collision vertex are 
retained in the analysis, thereby limiting the combinato- 
rial background from conversions and enabling a high 
detecting efficiency. The tracks are required to have at 
least 25 hits out of a maximum of 45 to avoid split tracks. 
Also, a minimum of 15 hits is required in the dE/dx mea- 
surement to obtain good dE/dx resolution [U [28j]. For 
particles directly originating from the collision, the col- 



lision vertex is added as an additional hit to further im- 
prove the momentum measurement [28j . 

Figure [1] panel (a) shows the 1//3 from TOF of parti- 
cles versus momentum in p + p collisions while panels (b) 
and (c) show the normalized dE/dx (na e ) distribution 
from the TPC as a function of px, without and with a 
requirement of high velocity 1 1//3 — 1 1 < 0.03, respectively. 
The quantity na e is defined as: na e — \n(dE / dx/ I e ) / R e , 
where dE/dx is the measured specific energy loss of a 
particle, and I e is the expected dE/dx of an electron. 
R e is the \n(dE/dx/I e ) resolution of an electron and is 
better than 8%. Electron candidates whose na e falls be- 
tween the lines indicated in Figure [T] panel (c) are re- 
tained in this analysis. With a perfect calibration, the 
na e for single electrons should follow a standard normal 
distribution. Figure [5] shows the na e distribution for 
0.4 <p T < 0.5 GeV/c after the cut of |l//3 - 1| < 0.03. 
The two dashed lines perpendicular to the x-axis repre- 
sent the range of the na e cut in this px region. A Gaus- 
sian plus exponential function, respectively representing 
the electron and hadron components, is used to fit the 
na e distribution. From the fit, we derive the purity and 
the na e cut efficiency on electron candidates as a function 
of pt, as shown in Fig. [3J The purity is defined within 
a range of na e (i.e. between the vertical dashed lines 
in Fig. [2]) as being the ratio of the electron counts in the 
area of the dashed Gaussian to the counts of all particles. 
The efficiency is defined to be the ratio of the electron 
counts under the dashed Gaussian within a range of na e 
to the total electron counts under the dashed Gaussian. 
The errors on the efficiency and purity are determined by 
adjusting the fit range. The electron yields are sensitive 
to the fit range since the hadron contamination increases 
in the smaller na e region thereby leading to large errors 
for the efficiency for pt > 0.8 GeV/c. Our exponential ex- 
trapolation of the na e distribution for the hadron compo- 
nent tends to over-estimate the background and therefore 
should be taken as an upper limit on the hadron contam- 
ination. By combining the velocity (/3) information from 
the TOF and the dE/dx from the TPC, electrons can 
be clearly identified from low to intermediate px (pt < 3 
GeV/c) for |77|<1 (H. 



B. Di-electron invariant mass distribution and 
background subtraction 

With a high purity for the electron samples, the e + 
and e~ from the same events are combined to gener- 
ate the invariant mass distributions (M ee ) of e + e~ pairs 
called unlike-sign distributions. The unlike-sign distribu- 
tions contain both signal and background. The signals 
are di-electrons from light- meson decays and heavy-flavor 
decays (charmonium and open charm correlation). The 
background results from random combinatorial pairs and 
correlated pairs. Electron candidates are required to be 
in the range of \rj\ < 1 and px > 0.2 GeV/c while e + e~ 
pairs are required to be in the rapidity range of \y ee \ < 1- 
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FIG. 1: (Color online) Panel (a): 1/(3 versus momentum of tracks from the TOF with \r)\ < 1 from 200 GeV p + p collisions. 
The line indicates the cut of — 1| <0.03. Panel (b): The normalized dE/dx distribution from the TPC as a function of pr- 
Panel (c): The normalized dE/dx distribution from the TPC as a function of pr with the cut of — 1| <0.03. An electron 
band is prominent, indicated by the lines, with the requirement of velocity close to the speed of light from TOF measurement. 




FIG. 2: (Color online) The na e distribution in 0.4<p T <0.5 
GeV/c after the cut of - 1| < 0.03 is applied. The solid 
curve represents a Gaussian plus exponential fit to the na e 
distribution. The dot-dashed line is for the hadron compo- 
nent and the dashed is for the electron contribution. The two 
dashed lines perpendicular to the x-axis represent the range 
of the n<r e cut to ensure a high purity for electron candidates 
in this pr range. 



The following two techniques are used for background 
estimation. In the like-sign technique, electron pairs 
with the same charge sign are combined from the same 
events. In the mixed-event technique, unlike-sign pairs 
are formed from different events. In order to ensure the 
events used in mixing have similar geometric acceptance 
in the detector, we only mix events which have collision 
vertices within 5 cm of each other along the beam line 
direction. 

Neither method represents the background perfectly. 
In the low-mass region, there is a correlated cross pair 
background (coming from two pairs of e + e~ from the 
same meson decays: Dalitz decays followed by a conver- 
sion of the decay photon, or conversions of multiple pho- 
tons from the same meson). This background is present 
in the like-sign distribution but not in the mixed-event 
background. On the other hand, due to the sector struc- 
ture of detectors and the different curvatures of positively 
and negatively charged particle tracks in the plane per- 
pendicular to the magnetic field, like-sign and unlike-sign 
pairs will have different acceptance. Moreover, in the 
high invariant mass range, there may be contributions 
from jet correlations which are absent from the mixed- 
event technique. 

Figure 21 shows the invariant mass distribution for 
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FIG. 3: (Color online) The purity and the no e cut efficiency 
for electron candidates as a function of pr in ] -7-7 1 < 1 for p + p 
collisions at y/s = 200 GeV. The squares represent the pu- 
rity and closed circles represent the na e efficiency. The pr 
positions of the last three data points for the efficiency are 
slightly shifted for clarity. The error bars are the quadrature 
sum of statistical and systematic uncertainties. A function 
form of A + B/[C + exp(D x pr)] is used to fit the efficiency 
data points and for the efficiency correction, in which A, B, 
C and D are the fit parameters. A constant component in 
the fit for pr> 0.8 GeV/c is driven by the study in Ref. [23| . 
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FIG. 4: (Color online) The electron-pair invariant mass distri- 
butions for unlike-sign pairs, like-sign, and mixed-event (com- 
binatorial) background in minimum-bias p + p collisions. The 
electron candidates are required to be in the range of < 1 
and have a pr greater than 0.2 GeV/c. The ee pairs were 
required to be in the rapidity range of \y e e\ < 1. The un- 
even bin widths are used based on the yields and the signal 
to background ratios. 



unlike-sign pairs, like-sign pairs, and mixed-event back- 
ground. The mixed-event distribution is normalized by 
a constant to match the like-sign distribution in the 
mass range 0.4-1.5 GeV/c 2 . In our analysis we subtract 
the like-sign background for M ee < 0.4 GeV/c 2 and the 
mixed-event background in the higher-mass region. For 
a cross-check on the consistency of the two methods we 
compare their shapes in the higher-mass region (Fig.[S]). 
A constant fits the ratio of like-sign over mixed-event dis- 
tributions for M ee >0.4 GeV/c 2 with X 2 /NDF = 14/15, 
as shown in the ratio plot in Fig. [5] (a). We also find 
that for electron pairs with significantly higher trans- 
verse momentum (as determined from the Barrel Electro- 
Magnetic Calorimeter triggered events), the shapes of 
like-sign and mixed-event distributions agree in the mass 
region of 1-3 GeV/c 2 [Ij. Fig. (b) shows the distri- 
bution of the difference of the azimuthal angles (A0) of 
the two electrons in the unlike-sign, like-sign and mixed- 
event pairs. The difference between like-sign and mixed- 
event pairs at low A</> can be attributed to cross-pair 
contributions. For M ee > 0.4 GeV/c 2 , the A0 distri- 
butions of like-sign and mixed-event pairs match nicely 
with each other as shown in Fig. [5] (c), indicating that 
mixed-event background subtraction is valid in the cor- 
responding mass region. 

As an additional check with a different method, we 
also perform the analysis by subtracting the like-sign 
background in the whole mass region. The difference 
in the di-electron yields from the mixed-event and like- 
sign methods is found to be within errors for M ee > 0.4 
GeV/c 2 . We also correct for acceptance differences be- 
tween the like-sign and unlike-sign electron pairs in both 
methods and will discuss the details in Sec. IIV CI 

Figure [5] (a) shows the di-electron continuum after 
background subtraction from both like-sign and mixed- 
event methods in p + p collisions at \/s = 200 GeV. The 
measurements are done requiring |y e + e -| < 1, \r] e \ < 1 
and pr(e) > 0.2 GeV/c and no efficiency correction has 
been applied. The two methods give consistent results 
for M ee > 0.4 GeV/c 2 . For the following, we use the 
yields from like-sign method for M ee < 0.4 GeV/c 2 and 
results obtained from the mixed-event method at higher 
mass as the default since the mixed-event background 
distribution matches the like-sign distribution and has 
better precision for M ee > 0.4 GeV/c 2 . The signal to 
background ratio in p + p collisions versus M ee is shown 
in Fig. El (b). 



IV. DI-ELECTRON CONTINUUM IN STAR 
ACCEPTANCE 

A. Cocktail simulation 

The di-electron pairs may come from decays of the 
light-flavor and heavy-flavor hadrons. They include 7r°, 
t], and rj' Dalitz decays: ir° —> r ye + e", r\ 



"fe + e 



and rj' -+ r ye + e ; vector meson decays: u — > 7r°e + e 
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FIG. 5: (Color online) Panel (a): The ratio of like-sign to mixed-event distributions in minimum- bias p + p collisions. Panel (b): 
The distribution of the difference of the azimuthal angles (A(f>) of the two electrons in the unlike-sign, like-sign and mixed-event 
pairs in minimum-bias p + p collisions. Panel (c): The A(f> distributions of unlike-sign, like-sign and mixed-event pairs for 
M ee >0.4 GeV/c 2 in minimum-bias p + p collisions. Errors are statistical. 




FIG. 6: (Color online) Panel (a): The di-electron (e + e~) continuum after background subtraction without efficiency correction 
in yfs = 200 GeV minimum-bias p + p collisions. Two methods of obtaining the background are indicated. Errors are statistical 
only. Panel (b): The signal over background ratio, plotted as a function of Af ee for NSD p + p collisions at \fs = 200 GeV. 
Errors are statistical. 



u! — » e + e~, p° — > e + e~, (j) ~ * r/ e+e 7 4> ~~ * e + e~, and 
J/ip — » e + e~; heavy-flavor hadron semi-leptonic decays: 
cc — ¥ e + e~ and bb -+ e + e~; and Drell-Yan contributions. 
We fit the invariant yields of mesons, previously mea- 
sured at STAR and PHENIX as discussed below, with 
the Tsallis functions (3l[, as shown in Fig. [7] (a). We use 
the Tsallis functions as input to a detector simulation 
in which the particles are decayed into di-electrons with 
the appropriate B.R.s. This GEANT detector simula- 
tion [32| uses a detailed description of the STAR geome- 
try in 2009. Simulated e + e _ cocktails from the various 
contributing sources, as shown in Fig. [7](b), are selected 
using the same cut conditions as those in the analyses 
of real events. The di-electron contributions from the 
gamma conversion 7 — > e + e~ in the detector material 
are accepted in both data and simulation subject to the 
same analysis cuts as well. The imperfect description 



of the material in the simulation leads to 3% system- 
atic uncertainty for the cocktail simulation for M ee <0.1 
GeV/c 2 . The Dalitz decays of?/ -> j°e + e~, uj -> 7r°e+e~ 
and 77' — > 7e + e~ are obtained using the Kroll-Wada ex- 
pression |33j |. 

For the Dalitz decays of 7r°, 77 and 77', we use the for- 
mula 

dN I imlf 2ml\ 1 / M e 2 e \ 3 2 a 

w 

in which, m e is the electron mass, M ee is the di-electron 
mass, and M% is the mass of the hadron which decays into 
the di-electron. F(M^ e ) is the electro- magnetic transition 
form factor. For the Dalitz decays of vector mesons cj and 
4> {A — > Be + e~), the formula is 
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FIG. 7: (Color online) Panel (a): The invariant yields of measured mesons fit with the Tsallis functions [3lJ in p + p collisions 
at yfs — 200 GeV. See text for details. Panel (b): The simulated raw di-electron continuum within STAR acceptance for 
y^s = 200 GeV minimum-bias p + p collisions. Different cocktail contributions are shown. 
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in which Ma and Mb are the mass of particle A and B, 
respectively. For all Dalitz decays except 77', the form 
factor is parameterized as 



F(Ml) 



1 



1 - MIA- 2 ' 



(3) 



in which A 2 is the form factor slope. For 77', we use the 
parametrization from 17] : 



\F(M 2 e 



2 \|2 



(l-M 2 e A- 2 ) 2 +T 2 A- 



(4) 



where the T 2 , is 1.99 x 10" 2 (GeV/c 2 ) 2 . The A parameters 
are listed in Table HI 

The p° — > e + e~ line shape is convoluted with the 
Boltzmann phase space factor [§3, HH and given by: 



dN 



dM ee d PT (m 2 - M 2 e ) 2 + m 2 (r^ + r ee r 2 )< 



_ M p fM 2 „ -AM, 
1 r 



2\ 3/2 



M ee V M 2 - AM 2 



x PS, 
(5) 

(6) 



r M p 

1 1 



(M 2 



Ami 



M ee V M p - 4m e 



1/2 



PS = 



M,. 



VM 2 e +p 2 T 



exp 



VmI+p't 



T 



(7) 



(8) 



where M p is 776 MeV/c 2 , M„ is the 7T mass, To is 149 
MeV/c 2 , T 2 is the B.R. of p° -+ e+e~, PS is the Boltz- 
mann phase space factor, and the inverse slope parameter 
T is 160 MeV [33 ]. We neglect any contribution from the 
interference effect between p° and u>. 

The invariant yield of it is taken as the avera ge o f 7r + 
and 7T" [H, [il|. The yields of 6 [H and are 
from STAR while the 77 [11], w [3] and J/t/> [451] yields 
are measurements by PHENIX. Table U lists the total 
yields at mid-rapidity {dN/dy\ y=0 ) in 200 GeV NSD p+p 
collisions. 

The last input we consider in the cocktail simulation 
is the cc cross section, which has been constrained by 
the published low-p^ D° spectrum in 200 GeV d+Au 
collisions . the non-photonic electron spectrum in 
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TABLE I: The total yields at mid-rapidity (dN/dy) from the Tsallis fit, decay B.R.s, and A parameters of hadrons in NSD 
p + p collisions at y/s = 200 GeV. 



meson 


dN 
dy 


relative uncertainty 


decay channel 


B.R. 


A (GeV/c ) 




1.28 




14% 


■ye + e~ 


1.174 x 10"^ 


1.756 ± 0.022 [36] 


V 


1.7 x 10" 


i 


.^o /o 


76 & 


7 n v i n — 3 


1 4- fl 1 8 \17] 
i.yo u.io 


p 


2.2 x 10" 


i 


15% 


e + e~ 


4.72 x 10~ 5 




UJ 


1.3 x 10" 


i 


21% 


e + e~ 


7.28 x 10" 5 




UJ 








ir°e+e- 


7.7 x 10~ 4 


2.24 ± 0.06 [37] 


<t> 


1.7 x 10" 


2 


20% 


e + e~ 


2.954 x 10~ 4 




<t> 








■qe + e~ 


1.15 x 10~ 4 


3.8 ± 1.8 [38] 


rf 


4.1 x 10" 


2 


29% 


7e + e~ 


4.7 x 10 -4 [39J 


1.8 ±0.4 [17] 


J/i> 


2.4 x 10" 


5 


15% 


e + e~ 


5.94 x 10" 2 





200 GeV p + p collisions [46j] , and the di-electron con- 
tinuum in this analysis. The details of these constraints 
will be shown in Sec. IIV CI The e + e~ shapes from open 
heavy-flavor pairs are obtained using PYTHIA6.416, in 
which the k T factor is set by PARP(91)=1 GeV/c, and 
the parton shower is set by PARP(67)=1 [47j]. With these 
parameters chosen, PYTHIA can describe the shape of 
STAR measured D° spectrum and non-photonic electron 
spectrum. The total contribution from the simulation 
is shown as the black solid curve in Fig. [7] (b). In the 
intermediate mass region, the di-electron continuum is 
dominated by the cc contribution. 



B. Efficiency and acceptance correction 

For the di-electron continuum, the efficiency cor- 
rections are applied within the STAR acceptance of 
\y e + e - \ < 1) \Ve\ < 1 and pr(e) > 0.2 GeV/c. The sin- 
gle electron efficiency includes the TPC tracking effi- 
ciency, TOF acceptance and detector response, and the 
dE/dx efficiency. Single electron tracking efficiency and 
acceptance in the TPC are determined by Monte Carlo 
GEANT simulations. The TOF acceptance and response 
efficiency for electrons is found to be 46% independent of 
Pt for 1 77 1 < 1 [25j . The efficiency of the na e cut, used to 
ensure a high purity for the electron sample, is close to 
100% at low p T and falls to -80% by p T = 0.8 GeV/c, 
as shown in Fig. [3J 

Figure [5] shows the efficiency for single electrons in the 
pseudo-rapidity range of < 1 in p + p collisions at 
y/s = 200 GeV. Open circles represent the TPC track- 
ing efficiency alone. Including TOF matching and re- 
sponse decreases the efficiency as shown by the triangles. 
With additional dE/dx cuts, the final efficiency for sin- 
gle electrons for |?y e | < 1 is shown as squares. For this 
analysis, 72% out of a possible 120 total TOF trays were 
installed and the efficiency dependence on azimuthal an- 
gle is shown in Fig. [5] for positive and negative r\ regions. 
We have accounted for the incomplete TOF acceptance 
in determining efficiencies for the di-electron spectra. 

The efficiency factor for the di-electron continuum 
within STAR's acceptance is obtained in two steps. We 
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FIG. 8: (Color online) The efficiency for single electrons as 
a function of pr in the pseudo-rapidity range of \r/\ < 1 from 
p + p collisions at y/s — 200 GeV. Open circles represent the 
TPC tracking efficiency. With additional TOF matching and 
response, the efficiency is shown as triangles. With additional 
dE/dx cuts, the final efficiency is shown as squares. 
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FIG. 9: (Color online) The TOF efficiency as a function of 
azimuthal angle in the positive and negative r\ region. The 
hole in the negative r\ region corresponds to the trays (28%) 
missing from the 2009 run installation. The 3 trays around 
270 degree in the positive r\ region have smaller matching 
efficiencies but have similar performances in other aspects. 
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obtain the input cocktail A within STAR acceptance by 
the method described in Sec. IIV Al The input cocktail 
includes the radiation energy loss and momentum reso- 
lution determined from GEANT simulations. The result 
is shown by the solid line in Fig. [10] (a). Then we ob- 
tain cocktail B from GEANT simulations with proper 
efficiency factors including the TPC tracking, TOF ac- 
ceptance and response, and dE/dx cut for single electrons 
as described above, and shown as the dashed line in Fig- 
ure [TU] (a). The ratio of these two is the efficiency factor 
for the di-electrons shown in Fig.[TU](b) ior p+p collisions 
at y/s = 200 GeV. The uncertainty on the efficiency fac- 
tor is about 10% with a negligible pr dependence. The 
final di-electron continuum within the STAR acceptance 
is obtained after this correction is applied and is discussed 
in Sec. HVCl 



C. Results 

The systematic uncertainties on the di-electron contin- 
uum are dominated by background subtraction (accep- 
tance difference between like-sign and unlike-sign elec- 
tron pairs and normalization of mixed-event distribu- 
tions) and hadron contamination. The acceptances of the 
like-sign and unlike-sign distributions are the same within 
5% for M ee > 0.1 GeV/c 2 due to the azimuthal symme- 
try and the solenoidal magnetic field. The small accep- 
tance differences due to track merging, sector boundaries, 
and dead channels have been evaluated using the differ- 
ences of the unlike-sign and like-sign distributions from 
the mixed-event technique. This difference is included in 
the systematic uncertainty. In addition, for M ee > 0.4 
GeV/c 2 , the normalization factor between the mixed- 
event and the like-sign distributions contributes 0-7% to 
the overall uncertainty. The systematic uncertainty from 
efficiency factors is about 10% with no significant mass 
dependence. The uncertainties in hadron contamination 
(hadrons from resonance decays mis-identified as elec- 
trons) are 0-32% and are mass dependent. Figure [TT] 
shows the relative systematic uncertainties from differ- 
ent sources for each mass bin. The normalization un- 
certainty in NSD p + p collisions is 8% [46[ . Additional 
contribution from the normalization uncertainty in VPD 
triggered minimum-bias events taking into account the 
trigger bias and vertex finding efficiency is 8% as deter- 
mined from PYTHIA simulations. The total normaliza- 
tion uncertainty for di-electron mass spectra is 11% in 
p + p collisions. Table [TT1 shows the systematic uncertain- 
ties from different sources. 

The uncertainties on the cocktail simulations include 
the decay form factors and the measured cross section for 
each hadron. By fitting the di-electron continuum, open 
charm [2^] and non-photonic electron spectra 46 1 simul- 
taneously, we find that the cc cross section in 200 GeV 
p + p collisions is 0.92 ± 0.10 ± 0.26 mb, consistent with 
earlier RHIC measurements [HI, HH . This is used to gen- 
erate the charm component in this paper. The systematic 



TABLE II: Systematic uncertainties from different sources for 
di-electron continuum. 



source 


contribution factors 


background subtraction 


0-27% 


contamination 


0-32% 


efficiency 


10% 


total normalization 


11% 


cocktail simulation 


14-33% 



uncertainties for the cocktail simulation as a function of 
mass are shown in Fig. 111! Future precise measurements 
of di-electron continuum in the intermediate mass region 
can further constrain the charm production mechanism 
in p + p collisions. 

After the efficiency correction, the di-electron contin- 
uum within the STAR acceptance is shown in Fig. [12] for 
= 200 GeV NSD p+p collisions. The di-electron mass 
spectrum is not corrected for momentum resolution and 
radiation energy loss effects. The ratio of data to cock- 
tail simulation is shown in the lower panel of Fig. 1121 
Within the uncertainties, the cocktail simulation is con- 
sistent with the measured di-electron continuum. The 
X 2 /NDF between data and cocktail simulation are 21/26 
for M ee >0.1 GeV/c 2 and 8/7 for l.K M ee < 3.0 GeV/c 2 . 
In the mass region 0.2 < M ee < 0.8 GeV/c 2 , the cocktail 
simulation is systematically higher than the measured 
di-electron continuum. However, they are also consistent 
within uncertainties. We find that better agreement be- 
tween the cocktail simulation and data can be achieved 
by applying an additional scale factor (56%) to the r] 
Dalitz decay contribution. Further details on this decay 
can be found in Sec. IVII 



V. VECTOR MESON PRODUCTION 

The yields of the w, 4> and J/ip long-lived vector mesons 
can be extracted from the di-electron continuum. We 
use the mixed-event technique to reconstruct the combi- 
natorial background beneath the respective peaks. The 
mixed-event distribution is normalized by a constant to 
match the like-sign distribution in the mass region of 0.4- 
1.5 GeV/c 2 , as discussed in Sec. lIIIBl The background is 
then subtracted to obtain the signal, which will still con- 
tain some residual background of di-electron pairs from 
other sources as described in Sec. IIV Al 

A two-component fit is used to extract the raw signal 
lj — > e + e~ from the residual background in the invari- 
ant mass range of 0.7 < M ee < 0.85 GeV/c 2 . The first 
component represents the line shape (the to — > e + e~ sig- 
nal shape), the second the residual background. The line 
shape of the w — > e + e~ invariant mass distribution, and 
the shape and magnitude of the background are deter- 
mined from the simulations described in Sec. IIV Al The 
systematic uncertainties of the w — > e + e~ raw yields are 
derived by changing the magnitude of the background al- 
lowed by the uncertainties of the cocktail simulation and 
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FIG. 10: (Color online) (a) The M ee distribution within STAR's acceptance {\y e + e - \ < 1, |ry e | < 1, and pr(e) > 0.2 GeV/c) 
from simulation. The solid line represents the input cocktail. The dashed line represents the cocktail from GEANT simulations 
taking into account the proper efficiency factors described in the text, (b) The efficiency of the di-electron continuum as a 
function of M ee within the STAR acceptance in p + p collisions at yfs = 200 GeV. 
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FIG. 11: (Color online) Systematic uncertainties as a function 
of mass from different source contributions. Also shown are 
the systematic uncertainties of the cocktail simulation. 
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the analysis cuts. The total contribution to the raw yield 
is about 20%. Figure [T3l shows the fit to the M ee distri- 
bution for 0.7 <M ee < 0.85 GeV/c 2 in three different p T 
bins for p + p collisions at y/s — 200 GeV. 

In order to present the final differential invariant cross 
section as a function of pr, the raw vector meson yield 
(uj — > e + e~) is corrected for acceptance and all detec- 
tor effects which reduce the measured raw yield relative 
to the actual yield. The total efficiency correction for 
w — > e + e~ for \y\ < 1 is shown in Fig. [14] We use the 
simulations described in Sec. IIV Al to determine this cor- 
rection, which accounts for limits in TPC acceptance and 
inefficiencies in TPC tracking, limits in the TOF accep- 
tance/response, and the rejection of signal due to the 
dE/dx cut. The invariant yield is defined as follows: 

d 2 N = N raw x NORM 

2npTdpTdy 2-rtpTdpTdy x N event x e x B.R. ' 

where N raw represents the raw signal counts, N event is 



FIG. 12: (Color online) The comparison of the di-electron 
continuum between data and simulation after efficiency cor- 
rection within the STAR acceptance in y/s = 200 GeV NSD 
p + p collisions. The di-electron continuum from simulations 
with different source contributions is also shown. The statisti- 
cal errors on the data are shown as bars, while the systematic 
uncertainties are shown as boxes. The 11% normalization 
uncertainty on the data is not shown. The band on top of 
the solid curve illustrates the systematic uncertainties on the 
cocktail simulation. 
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FIG. 13: (Color online) The M ee distribution in the range 0.65 < Mee < 0.95 GeV/c 2 for three different Pt bins after the 
mixed-event background subtraction in p + p collisions at y/s = 200 GeV. The curves represent fits in the range 0.7 < Mee < 0.85 
GeV/c 2 as described in the text. Errors on data points are statistical. 
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FIG. 14: (Color online) The efficiency including STAR accep- 
tance for lo — > e + e~ for \y\ < 1 in p + p collisions at yfs — 200 
GeV. 



the total event number, e is the total efficiency and ac- 
ceptance correction factor, B.R. is the branching ratio 
for to — > e + e~, and NORM is the normalization factor 
(64 ± 5% for VPD triggered minimum-bias events taking 
into account the trigger bias and vertex finding efficiency, 
determined by, and corrected for, using PYTHIA simu- 
lations) . The oj invariant yield in NSD p + p collisions at 
y/s = 200 GeV is presented in Fig. [15J The systematic 
uncertainties are dominated by uncertainties in the two- 
component fit, and uncertainties in total efficiency which 
are also described in Sec. I IV Al Table Qll] lists the de- 
tailed systematic uncertainties for lo yields from different 
sources. Our lo yields from di-electron decays are consis- 
tent with previous results 49] and with a prediction from 
a Tsallis fit, which fits spectra of other particles and high 



TABLE III: Systematic uncertainties from different sources 
for lo yields. 



source 


contribution factors 


two component fit 


9-16% 


na e cut 


3-7% 


efficiency 


10% 


total normalization 


11% 



Pt uj yields [U . We obtain a mid-rapidity dN/ dy of 
0.10 ± 0.02 (stat.) ± 0.02 (sys.) for the lo. 

In addition to lo — > e + e~ yields, we also obtain the 
mid-rapidity yields, dN/dy, for the <j> and J/tp particles 
in NSD p + p collisions at y/s = 200 GeV. Due to limited 
statistics for each particle, the invariant mass signal can 
only be extracted over allpr, rather than individual bins. 
As before, we use a two-component fit in the region of 
M ee for 0.98 < M ee < 1-04 GeV/c 2 for the -> e+e" , and 
3<M ee <3.16 GeV/c 2 for J/i/> -> e+e~. The line shapes 
of <j) — > e + e~ and J/ip — > e + e~ arc from simulations as 
discussed in Sec. IIV Al The residual background shape 
and magnitude are obtained from the cocktail simula- 
tion, shown in Sec. IIV Al The systematic uncertainties of 
the <j) and J/tp dN/dy due to the residual background are 
evaluated by changing the magnitude of the background 
allowed by the uncertainties of the cocktail simulation. 
The total efficiency correction for each particle is evalu- 
ated in the same way as for the lo — > e + e~ analysis. Since 
the correction depends on px, we calculate a weighted 
average using the predicted spectra from the previously 
mentioned Tsallis fit as px weights. The total efficiency 
and acceptance correction factors for <j> and J/ifj are 4.4% 
and 3.4%, respectively. The systematic uncertainty on 
the total efficiency correction is estimated to be 10%. 
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FIG. 15: (Color online) The uj — > e + e" invariant yield, di- 
vided by its B.R., in NSD p + p collisions at s/s = 200 GeV. 
The open circles represent PHENIX published results [49l |. 
The bars are statistical errors and boxes are systematic un- 
certainties. A normalization uncertainty of 11% is not shown. 
The yields at the center of the pr bin are corrected for finite 
bin width. The line represents the yields of lj from Tsallis 
function fit to high pr CJ yields measured from its hadronic 
decay and described in Sec. IIV Al 



TABLE IV: Systematic uncertainties from different sources 
for <f> and J/if) dN/dy. 



source 


contribution factors for <f> 


for J /if) 


two component fit 


10% 


7% 


na e cut 


8% 


27% 


efficiency 


10% 


10% 


total normalization 


11% 


11% 



The normalization uncertainty is 11% in p + p collisions. 
Table IIVI lists the detailed systematic uncertainties for 
the (f> and J/ip dN/dy from different sources. Figure [TBI 
shows the fits to the M ee distributions used to obtain the 
mid-rapidity yields dN/ dy for <fi and J/tf> in NSD p + p 
collisions at y/s — 200 GeV. The final yields in those fits 
are subject to the total efficiency correction. The dN/dy 
for <j> is 0.010 ± 0.002 (stat.) ± 0.002 (systL consistent 
with the measurements from 4> — > K + K~ [23, The 
dN/dy for J/tp is (2.1 ± 0.7 (stat.) ± 0.7 (systV) x 
10~ 5 , consistent with previous measurements [45l. l51j. 



VI. THE RARE DECAY: r\ e+e" 

As discussed in previous sections and shown in Fig. [12] 
the cocktail can describe the data reasonably well around 
the 77 mass without the r/ — > e + e~ decay channel. We 
zoom into the low-mass range and show the data and 
cocktail comparison in Fig. [T7J The dot-dashed peak is 
the 1] — > e + e~ contribution with the upper limit of its 
B.R. from the PDG [H, which is 2.7 x 10~ 5 [52]. The 
dashed curve is a two-component fit with the 77 — » e + e~ 
decay channel included. The fit function is 

AxjV, + JJx Cocktail, (10) 

in which N v is the expected 77 contribution with the line 
shape of 77 — > e + e~ after detector simulation, Cocktail is 
the expected cocktail contribution described in Sec. IIV Al 
without 77 — > e + e~, and A and B are fit parameters. The 
A and B represent the B.R. for the 77 — > e + e~ and a scale 
factor for the cocktail contribution, respectively. The 
solid curve shown in Fig. ll7l represents the 77 — > e + e~ con- 
tribution from the fit. It gives the B.R. of the 77 — > e + e~ 
to be (-9.6 ± 5.9 (stat.) ± 5.3 (syst.)) x 10~ 6 . The 
negative value of B.R. of the 77 — > e + e~ is due to the 
statistical fluctuation at M ee — 0.55 GeV/c 2 . The sys- 
tematic uncertainties are dominated by background sub- 
traction (34%, the difference between mixed event back- 
ground subtraction and like-sign subtraction), electron 
purity (31%, determined by varying the no~ e cuts), and 
track quality cut (27%, determined by changing the cut 
of distance of closest approach between the track and the 
collision vertex). 

In addition, although the Dalitz decay yield of ?; — > 
7e + e~ is consistent with the cocktail expectation from 
the Tsallis fit to rj measurements for pt > 2 GeV/c de- 
scribed in Sec. IIV Al the nominal value from the fit in 
Fig.[l7]is about 56% of the input cocktail (B=56%). This 
additional factor leads to a better agreement between 
data and cocktail simulation compared to that shown in 
Fig. [T3J This is equivalent to a smaller value of N v and 
has to be taken into account when the upper limit at the 
90% Confidence Level (CL) for the B.R. of 77 ->• e+e" is 
estimated. 

With different background subtraction, electronic pu- 
rity, and track quality cuts, we repeat the fit process de- 
scribed above to obtain the parameter A, the B.R. of the 
77 — > e + e~ . The difference of the A values is attributed 
to the systematic uncertainties. In this fit procedure, we 
find that cut conditions do not contribute to the point-to- 
point variation around 77 mass range and the statistical 
error of the parameter A remains unchanged. Therefore, 
the significance of an observable signal above the back- 
ground only depends on the statistical fluctuation. To 
estimate the upper limit for the B.R. of 77 — > e + e~ at 
the 90% CL due to a possible statistical fluctuation, wc 
utilize the statistical error 5.9 x 10 -6 for the B.R. of the 
77 — > e + e _ from the fit and set the lower bound of the 
B.R. of 77 — > e + e~ to be zero instead of the negative value 
from the fit. The upper limit is found to be 5.9 x 10~ 6 
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FIG. 16: (Color online) Panel (a): The M ee distribution after the mixed-event background subtraction for 0.96 < M ee < 1.08 
GeV/c 2 in p + p collisions at yfs — 200 GeV. The curve represents the fit in the range of 0.98 < M ee < 1.04 GeV/c 2 . Errors on 
data points are statistical. Panel (b): The M ee distribution after the mixed-event background subtraction for 2.98 <M ee <3.22 
GeV/c 2 in p + p collisions at </s = 200 GeV. The curve represents the fit in the range of 3 < M ee < 3.16 GeV/c 2 . There is no 
count in the unlike-sign and like-sign electron-pair distributions for 3.02 < M ee < 3.04 GeV/c 2 due to statistical fluctuations. 
Errors on data points are statistical. 



x 1.65/0.56 = 1.7 x 10 -5 , in which 1.65 is the upper 
endpoint of a confidence interval (the lower endpoint is 
zero) with a 90% CL for a standard normal distribution. 



1,1, 

— T|H>ee 1 

- - - Cockts 
- — r|H>ee 


1 , 1 , 1 1 , | 1 , 

'DG upper limit 
il + T|— >ee fit 
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FIG. 17: (Color online) M ee distribution for di-electron pro- 
duction in p + p collisions at ^/s = 200 GeV. The dashed 
curve is the cocktail with the 77 — > e + e~ decay channel in- 
cluded in the fit. The 77 — > e + e~ contribution is shown as 
the solid curve. The dot-dashed peak is the 77 — > e + e~ con- 
tribution with the upper limit of its branching ratio from the 
PDG [H, which is 2.7 x 10~ 5 . 



These results provide a promising first glimpse of a pro- 
gram on searching for rare decays of hadrons produced 



in relativistic heavy-ion collisions at STAR. With large 
hadron yields, high efficiency for electrons at low momen- 
tum and high mass resolution, STAR provides a unique 
tool for such a program in the years to come. 

VII. SUMMARY 

The di-electron continuum is measured in i/s = 200 
GeV non-singly diffractive p + p collisions within the 
STAR acceptance. The cocktail simulations are consis- 
tent with the data and provide a reference for the future 
heavy-ion studies. The u) invariant yields are consistent 
with the previous publications. The dN/dy for <f> and 
J/ip are 0.010 ± 0.002 (stat.) ± 0.002 (syst.) and (2.1 ± 
0.7 (stat.) ± 0.7 (syst.)) x 10~ 5 , respectively. These re- 
sults are consistent with the previous measurements from 
4> — > K + K~ and J/ip — > e + e~. Our measurement low- 
ers the current world limit of the branching ratio of the 
n e+e- from 2.7 x 10~ 5 to 1.7 x 10~ 5 . 
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